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Abstract 

In low-density nuclear matter which is relevant to the crust region of neutron stars and 
collapsing stage of supernovae, non-uniform structures called "nuclear pasta" are expected. So 
far, most works on nuclear pasta have used the Wigner-Seitz cell approximation with anzats 
about the geometrical structures like droplet, rod, slab and so on. We perform fully three- 
dimensional calculation of non-uniform nuclear matter for some cases with fixed proton ratios 
and in beta-equilibrium based on the relativistic mean-field model and the Thomas-Fermi 
approximation. In our calculation typical pasta structures are observed. However, there 
appears some difference in the density region of each pasta structure. 

1 Introduction 

In 1934, W.Baade and F.Zwicky proposed the existence of neutron stars only one year after the 
discovery of neutron. Estimating the binding energy of neutron stars, they predicted that neutron 
stars are made by supernova explosions. Neutron star had been an imaginary object for 30 years, 
until pulsars were discovered by J. Belle and A.Hewish in 1967. Neutron star has a radius of about 
10 km and the mass about 1.4 times of the sun, and consists of four parts pQ. The region around 
0.3km from outside is called "outer crust", where Fe nuclei are expected to form a Coulomb lattice. 
Around 0.3 — 1km is called "inner crust" with a density about 0.3-0.5po- There are neutron-rich 
nuclei in lattice and dripped neutrons in a superfluid state. Two central regions with higher density 
are called "outer core" and "inner core" , where speculated are proton super-conductivity, neutron 
super-fluidity [2 , meson condensations [3 , hyperon mixture [3] , or quark states |5j [6] . 

The transition of matter composition with the change of density inside neutron stars causes 
a question: does it change continuously or suddenly? A sudden change of matter property is 
generally accompanied by a first-order phase transition which causes an appearance of the mixed 
phase. 

Ravenhall et al [7] suggested the existence of non-uniform structures of nuclear matter, i.e. the 
structured mixed phase. They suggested five types of structures as droplet, rod, slab, tube, and 
bubble. Due to its geometrical shapes which depend on the density, we call it "nuclear pasta" like 
spaghetti and lasagna etc [7J [8] . Many workers have suggested the existence of pasta structures 
in low-density nuclear matter, relevant to the crust region of neutron stars and the collapsing 
stage of supernovae. The existence of the pasta structures at the crust of neutron stars may not 
influence on the bulk property and structure of neutron stars. However, it should be important 
for the mechanism of glitch, the cooling process of neutron stars, and the thermal and mechanical 
properties of supernova matter. 

The pasta structures presented by Ravenhall have geometrical symmetries. So we can treat the 
system with the Wigner-Seitz (WS) cell approximation. Because of the convenience, the WS cell 
approximation has been very often used. But there may be some shortcomings. First, the existence 
of some kinds of structures other than the typical pasta structures were suggested. For example, in 
compressive process of matter, two droplets connect with each other and form dumbbell-like pieces 
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[9]. Other examples are double diamond [10j[TT] and gyroid [TT] shapes of matter suggested by using 
compressible liquid-drop model. These shapes can't exist as a ground state at zero temperature, 
so they might exist in supernova matter at finite temperatures. Such structures are impossible to 
be described by the WS approximation. It is worth trying to calculate without any approximation 
whether or not these structures exist as a ground state or an excited state. 



2 Method 

2.1 Relativist ic Mean Field Theory 

There are several approaches for studying pasta structures in the literature, as liquid-drop model [7, 
[TT] . Thomas- Fermi model [12J[T3], Hartree-Fock [14], quantum molecular dynamics (QMD) [T5] l9], 
relativistic mean field model (RMF) p~6j [J7] . In the studies using liquid-drop model and RMF 
model, always used was the Wigner-Seitz (WS) cell approximation, where only the typical pasta 
structures are considered. QMD calculation does not assume any specific structure of baryons. But 
uniform electron distribution, though it should be almost uniform, is assumed. The Thomas-Fermi 
calculation in Ref. [12] and HF calculation [14 used the periodic boundary condition and did not 
assume any geometrical symmetry in structure. However, the size of the periodic cell was not large 
enough for quantitative discussion. 

In this paper, we use the Thomas- Fermi approximation for baryons with interaction by the RMF 
model. This model is not only simple for numerical calculation but also quantitatively reliable for 
the properties of finite nuclei and the saturation property of matter. 

We start from the simple Lorentz-scalar Lagrangian with nucleons, electrons and a, a;, p- mesons 
as follows, 



L = i\) 



1 + 73 

n^d/i - m* N - guN^u^ - gpNj^r-R^ - e — - — j^A^ 
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where the effective mass is written as m* N = mjv — g a N^, and the potential energy of sigma meson 
U(a) = \bm N {g a af - \c{g a af. 

From the variational principle for this Lagrangian, we get equations of motion for nucleons, a, 
oj, p- mesons and the Coulomb potential as follows: 



1 + i~z 

i^df,, - m* N - guNl^Ufj, - gpN^r-R^ - e — - — 7*M M 



V> = o, 
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V 2 A = e 2 (fy°(l + r 3 )/2^ - &7o^e) = e 2 p ch . 



(4) 
(5) 
(6) 



By the mean-field approximation for meson fields and the static approximation for electric field, 
we have (a) = a, (oj^) = S^^ujq, (R%) = S^^Sa^Ro and V = (Aq). The bra-ket () represents the 
expectation value in the ground- state of nuclear matter, and we assume the space and time-reversal 
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2.1 Relativistic Mean Field Theory 



symmetries for this state, uj and p mesons can't have finite value in space components. Therefore 
we finally get the equations for mesons as follow. 



V 2 cr + m^(J 



dU 



da 



■9*b{W) 



(J) 

- V 2 uj + mluo = 9ujn{p p + pn) (8) 

- V 2 R + m 2 p R = g P N(p P - Pn) (9) 
By the Thomas-Fermi approximation at zero temperature, momentum distributions of nucleons 
have step-functional forms and plane- wave solutions for Dirac equation. Scalar density and density 
of nucleons p s = Ylt=i ^(^O^O^)? P — ^2f=i ^Pii 00 )^ &(%) are written using ^ (i = 1, 2, • • • A) as 
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,here we use E% = \/m^ + k 2 and the electron density is given as p e = —J d 3 x('0e7oV ; e) < 0. 
The energy density of nuclear matter is obtained by the energy-momentum tensor and putting 
(/i,*/) = (0,0), as 
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Therefore, we get the total energy 
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2.2 Numerical Calculation 



We use the same parameter set of Ref. [16] , which reproduces the saturation property of nuclear 
matter (Fig. [I]), and the binding energies and the proton mixing ratios of finite nuclei. We list the 
parameter set in Table HJ 



QaN guN QpN 


b 


c m a (MeV) 


rriu (MeV) 


m p (MeV) 


6.3935 8.7207 4.2696 


0.008659 


-0.002421 400 


783 


769 



Table 1: Parameter set 




Figure 1: (a) Binding energy and (b) pressure without the contribution of electrons. Red line 
corresponds to the case of proton mixing ratio Y p = 0.5 and green, blue, purple, sky blue and 
yellow correspond to 0.4, 0.3, 0.2, 0.1 and 0, respectively. 



2.2 Numerical Calculation 

To simulate infinite matter, we employ a periodic boundary condition to the calculation cell with 
a cubic shape. Desirable cell size is large as possible. We divide the cell into grid points. The 
density distributions and the meson-field profiles are represented by their local values on the grid 
points. Giving an average densities of protons, neutrons and electrons, but density distributions 
are randomly provided. The suitable density distributions and fields are searched for. 

The meson fields and the Coulomb potential are obtained by solving the equations of motion 
|3j)-(|6j). To obtain the density distributions of nucleons and electrons we introduce local chemical 
potentials. The equilibrium state is determined so that the chemical potentials are independent of 
the position. An exception is the region with a particle density equals to zero, where the chemical 
potential of that particle becomes higher. 

M P = \/ (3^p p (r)) 2/3 + (m* N (r)) 2 + 9uN lj (r) + g pN R (r) - V (r) (15) 

= \J (37r 2 p n (r)) 2/3 + (m* N (r)) 2 + g^N^o (r) - g P NRo (r) (16) 
Me = (3^ 2 Pe (r)) 1/3 + F(r) (17) 

Starting from given density distributions, we repeat the following procedures to attain unifor- 
mity of chemical potentials. 
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2.3 Coulomb Energy 



1. Calculate the chemical potentials on all of the grid points. 

2. Compare chemical potentials on the neighbor six grid points. 

3. If the chemical potential of the point under consideration 

is larger than that of another, give some part of the density to 
the other grid point. 



2.3 Coulomb Energy 

In the Coulomb and electron energies (fT4]h we consider the energy only within the calculating 
cell without interaction with the neighbors. Therefore we do not include the Coulomb energy 
of the higher order, although we solve the Poisson equation completely. The dipole interaction 
which occupies the biggest contribution among the higher order interactions, comes up to 5% of 
the Coulomb energy. So we subtract the dipole moment of the electric field by carrying out a 
translation of the coordinate so that the dipole moment of charge density in the cell diminishes. 

3 Results 

We calculate three-dimensional structures of low-density nuclear matter at zero temperature and 
obtain the energy or pressure as function of density, i.e. the equation of state (EOS). In this report, 
let us first demonstrate the cases with the fixed proton mixing ratio, Y p = N p /Nb, and then with 
beta-equilibrium. Particularly, we set the proton mixing ratio to Y p =0.5, 0.3 and 0.1. Because 
these cases should be the one of typical nuclear matter, relevant for the supernova core, and for 
protoneutron stars. 

In some of the previous works which have used the WS cell approximation, nuclear matter 
might be enforced to have some typical pasta structures. In reality, however, there might appear 
some intermediate shapes in the density regions where structure changes. With the WS cell 
approximation, the cell size which gives the minimum energy density has been calculated [16]. To 
save the computational effort, we make use of that size of the WS cell for the size of our cell with 
a periodic boundary condition. 




(a) droplet (b) rod (c) slab (d) tube (e) bubble (f) uniform 



Figure 2: Proton density distribution for symmetric nuclear matter (Y p = 0.5) 

Red color corresponds to the highest density 0.08 fm -3 , and blue corresponds to fm 
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Figure 3: EOS of matter with Y p = 0.5. (a) Energy per nucleon obtained by the present calculation 
(colored circles) compared with the results of Ref. [16 J. (b) Total pressure and (c) baryon partial 
pressure. Black dotted lines in (a), (b) and (c) show the cases of uniform matter. 



First, we show the result for symmetric nuclear matter with a proton mixing ratio Y p = 0.5 
in Fig. [2j Panels (a), (b), (c), (d), (e), and (f) correspond to droplet, rod, slab, tube, bubble and 
uniform, respectively. In our calculation, all of the typical pasta structures are seen. 

The binding energy, the total pressure and the baryon partial pressure are presented in Fig. [3l 
The line with colors, dashed line and dots correspond to the results with the WS cell approximation, 
the case of uniform matter, and our results by the three-dimensional calculation, respectively. The 
density region with numbers 1,2,- •• ,6 separated with vertical dashed lines indicate that the 
structures in Fig. [2] (a), (b), • • • , (f) appear. Note that the density range for each pasta structure 
is slightly different from the previous study with the WS cell approximation. 




(a) droplet 



(b) rod 



(c) slab 



(d) tube 



(e) uniform 



Figure 4: Same as Fig. [2] for Y p — 0.3. But the density range from blue to red is to 0.075 fm 
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Figure 5: Same as Fig. [3] for Y p = 0.3. 
Next, we show the result for proton mixing ratio Y p = 0.3. in Fig. [5j The cell size is set to be 
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the same as in Ref. [16]. The difference between two results is that the bubble structure does not 
appear in the present calculation. 

Comparing the cases of Y p = 0.5 and 0.3, the upper limit of the density where non-uniform 
structures appear is different. In the case of symmetric nuclear matter, the density region of non- 
uniform matter roughly corresponds to the spinodal region, where dP/dpB < 0, while the non- 
uniform region is slightly wider than the spinodal region for Y p = 0.3. This may be because the 
symmetric nuclear matter behaves congruently (phase transition in a single chemical component), 
while the liquid-gas mixed phase in matter with Y p = 0.3 is non-congruent and the values of Y p in 
two phases are generally different. Therefore the instability of uniform matter is determined not 
only by the compressibility of uniform matter but also by the chemical composition of the mixed 
phase to be formed. 

Next, we show the result of Y p = 0.1. 




(a) droplet (b) rod (c) slab 



Figure 6: Same as Fig. [2] for Y p = 0.1. The density range from blue to red is to 0.05 fm 




Baryon density [fm-3] Batyon Density [fm- 3 ] Baryon Density [fm" 3 ] 

(a) Binding energy (b) Total pressure (c) Baryon pressure 

Figure 7: Same as Figs. [3] and [5] for Y p = 0.1. 

In Fig. [7J we show the baryon density dependence of binding energy, total pressure, and baryon 
partial pressure. We use the lines and dots in Non-uniform structures of only droplet, rod, and 
slab appear. 
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(a) Y p = 0.5 



(b) Y p = 0.3 



(c) Y p = 0.1 



Figure 8: Comparison of the cell and the structure size with the WS approximation and three- 
dimensional calculation (Y p =0.3) 

Line with colors : WS approximation, Dots: Three-dimensional calculation. Upper line with color 
i?ceii : ce ll size, Bottom line with color i?d : Non- uniform structure size 

To compare quantitatively our calculation with the WS cell approximation, we show the cell 
and the structure size in Fig. [U Upper and lower lines represent the cell and the structure size in 
one cell, respectively. Non-uniform structure size is defined by way of a density fluctuations as 



We show in Fig. [8] the size of the cell and the non-uniform structure. The density where each 
pasta structure appears is mostly in agreement with Ref. [16]. However, there are large differences 
in the transient region of pasta structures. 

In our calculation, we set, as the initial condition, the density distributions of fermions (nucleons 
and electron) and mesons randomly. So, the converged density distribution sometimes are trapped 
in local minimum states particularly at the density region where structure changes. 

Let us discuss the case of beta-equilibrium at zero temperature, which is relevant to the realistic 
neutron star matter. 




, (for droplet, rod and slab) 
^^), (for tube and bubble) 



(18) 
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(a) WS approx. (b) proton (c) neutron (d) electron 



Figure 9: Density distributions of nucleons and electron in the case of beta-equilibrium. The left 
panel (a) shows the result reported in Ref. [16 with WS approximation. The horizontal axis is the 
distance from the center of the cell, and the vertical axis densities of nucleons and electron. 




' Baryon Density [fm-3] Baryon Density [fnv 3 ' ] Baryon Density [fm 3 1 

(a) Binding energy (b) Total pressure (c) Proton mixing ratio Y p 

Figure 10: Energy, pressure, and Y p of matter in beta-equilibrium. 

In Figs.[9]and[T0j we present density distributions by the WS cell approximation and the present 
three-dimensional calculation. While some pasta structures appeared for fixed proton mixing ratio, 
only the droplet structure appears for the beta-equilibrium case. This result is qualitatively the 
same as in the WS cell approximation. 

Nuclear matter in the crust region of neutron stars includes the very small fraction of protons 
and electrons. So the contribution to the energy and pressure comes almost from neutrons and a 
little from protons and electrons. Consequently the binding energy and pressure of uniform matter 
are almost the same as those of non-uniform matter. The proton mixing ratio Yp, however, is 
significantly different between uniform and no-uniform cases. Though it approaches to zero in the 
zero-density limit for uniform matter, there is a steep rise for non-uniform matter. It approaches 
to approximately 0.5, which corresponds to a neutral atom in the low-density limit. 

4 Conclusion 

We have performed the three-dimensional calculations for low-density nuclear matter and presented 
first results, based on relativistic mean-field theory and Thomas-Fermi approximation. For some 
fixed proton mixing ratios and for /^-equilibrium nuclear matter, we have obtained non-uniform 
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structures and the relevant EOS. The observed structures are typical "pasta" structures, which were 
previously studied with the WS cell approximation. For Y p — 0.5, all of the typical pasta structures 
appeared. For Y p = 0.3 and 0.1, however, some of the pasta structures did not appear. The density 
region for each type of pasta structure was also different between the WS cell approximation and 
the three-dimensional calculation. 

These may be due to the difference of the shapes of the cell. The WS cell for droplet, rod, slab, 
tube and bubble structures are sphere, cylinder, plate, cylinder and sphere, respectively. On the 
other hand, the cubic cell is used in the three-dimensional calculation. In any way, it is desirable 
to avoid the effects of the shape of the used cell. We therefore are planning to enlarge the cubic 
cell so that several periods of structures can emerge in it. 

For /^-equilibrium, we got only the droplet structure. This result is qualitatively the same as 
the one with the WS cell approximation. However, we cannot definitely state that nuclear matter 
in /^-equilibrium have only droplet, because other works with different interaction have observed 
pasta structures other than droplet. In /^-equilibrium, proton mixing ratio is largely different from 
that in uniform matter. It tends to zero for uniform matter at zero-density limit, while it increases 
in the case of non-uniform matter, and the asymptotic value should be 0.5 which represents atoms. 
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